Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases.

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Integrative analyses highlight functional regulatory variants associated with neuropsychiatric diseases.

Noncoding variants of presumed regulatory function contribute to the heritability of neuropsychiatric disease. A total of 2,221 noncoding variants connected to risk for ten neuropsychiatric disorders, including autism spectrum disorder, attention deficit hyperactivity disorder, bipolar disorder, borderline personality disorder, major depression, generalized anxiety disorder, panic disorder, post-traumatic stress disorder, obsessive-compulsive disorder and schizophrenia, were studied in developing human neural cells. Integrating epigenomic and transcriptomic data with massively parallel reporter assays identified differentially-active single-nucleotide variants (daSNVs) in specific neural cell types. Expression-gene mapping, network analyses and chromatin looping nominated candidate disease-relevant target genes modulated by these daSNVs. Follow-up integration of daSNV gene editing with clinical cohort analyses suggested that magnesium transport dysfunction may increase neuropsychiatric disease risk and indicated that common genetic pathomechanisms may mediate specific symptoms that are shared across multiple neuropsychiatric diseases.

Pan-conserved segment tags identify ultra-conserved sequences across assemblies in the human pangenome.

The human pangenome, a new reference sequence, addresses many limitations of the current GRCh38 reference. The first release is based on 94 high-quality haploid assemblies from individuals with diverse backgrounds. We employed a k-mer indexing strategy for comparative analysis across multiple assemblies, including the pangenome reference, GRCh38, and CHM13, a telomere-to-telomere reference assembly. Our k-mer indexing approach enabled us to identify a valuable collection of universally conserved sequences across all assemblies, referred to as "pan-conserved segment tags" (PSTs). By examining intervals between these segments, we discerned highly conserved genomic segments and those with structurally related polymorphisms. We found 60,764 polymorphic intervals with unique geo-ethnic features in the pangenome reference. In this study, we utilized ultra-conserved sequences (PSTs) to forge a link between human pangenome assemblies and reference genomes. This methodology enables the examination of any sequence of interest within the pangenome, using the reference genome as a comparative framework.

Mutations in human DNA methyltransferase DNMT1 induce specific genome-wide epigenomic and transcriptomic changes in neurodevelopment.

DNA methyltransferase type 1 (DNMT1) is a major enzyme involved in maintaining the methylation pattern after DNA replication. Mutations in DNMT1 have been associated with autosomal dominant cerebellar ataxia, deafness and narcolepsy (ADCA-DN). We used fibroblasts, induced pluripotent stem cells (iPSCs) and induced neurons (iNs) generated from patients with ADCA-DN and controls, to explore the epigenomic and transcriptomic effects of mutations in DNMT1. We show cell type-specific changes in gene expression and DNA methylation patterns. DNA methylation and gene expression changes were negatively correlated in iPSCs and iNs. In addition, we identified a group of genes associated with clinical phenotypes of ADCA-DN, including PDGFB and PRDM8 for cerebellar ataxia, psychosis and dementia and NR2F1 for deafness and optic atrophy. Furthermore, ZFP57, which is required to maintain gene imprinting through DNA methylation during early development, was hypomethylated in promoters and exhibited upregulated expression in patients with ADCA-DN in both iPSC and iNs. Our results provide insight into the functions of DNMT1 and the molecular changes associated with ADCA-DN, with potential implications for genes associated with related phenotypes.

Evaluation of C4 Gene Copy Number in Pediatric Acute Neuropsychiatric Syndrome.

Pediatric acute-onset neuropsychiatric syndrome (PANS) is an abrupt-onset neuropsychiatric disorder. PANS patients have an increased prevalence of comorbid autoimmune illness, most commonly arthritis. In addition, an estimated one-third of PANS patients present with low serum C4 protein, suggesting decreased production or increased consumption of C4 protein. To test the possibility that copy number (CN) variation contributes to risk of PANS illness, we compared mean total C4A and total C4B CN in ethnically matched subjects from PANS DNA samples and controls (192 cases and 182 controls). Longitudinal data from the Stanford PANS cohort (n = 121) were used to assess whether the time to juvenile idiopathic arthritis (JIA) or autoimmune disease (AI) onset was a function of total C4A or C4B CN. Lastly, we performed several hypothesis-generating analyses to explore the correlation between individual C4 gene variants, sex, specific genotypes, and age of PANS onset. Although the mean total C4A or C4B CN did not differ in PANS compared to controls, PANS patients with low C4B CN were at increased risk for subsequent JIA diagnosis (hazard ratio = 2.7, p value = 0.004). We also observed a possible increase in risk for AI in PANS patients and a possible correlation between lower C4B and PANS age of onset. An association between rheumatoid arthritis and low C4B CN has been reported previously. However, patients with PANS develop different types of JIA: enthesitis-related arthritis, spondyloarthritis, and psoriatic arthritis. This suggests that C4B plays a role that spans these arthritis types.

Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability.

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

Maternal attachment insecurity, maltreatment history, and depressive symptoms are associated with broad DNA methylation signatures in infants.

The early environment, including maternal characteristics, provides many cues to young organisms that shape their long-term physical and mental health. Identifying the earliest molecular events that precede observable developmental outcomes could help identify children in need of support prior to the onset of physical and mental health difficulties. In this study, we examined whether mothers' attachment insecurity, maltreatment history, and depressive symptoms were associated with alterations in DNA methylation patterns in their infants, and whether these correlates in the infant epigenome were associated with socioemotional and behavioral functioning in toddlerhood. We recruited 156 women oversampled for histories of depression, who completed psychiatric interviews and depression screening during pregnancy, then provided follow-up behavioral data on their children at 18 months. Buccal cell DNA was obtained from 32 of their infants for a large-scale analysis of methylation patterns across 5 × 106 individual CpG dinucleotides, using clustering-based significance criteria to control for multiple comparisons. We found that tens of thousands of individual infant CpGs were alternatively methylated in association with maternal attachment insecurity, maltreatment in childhood, and antenatal and postpartum depressive symptoms, including genes implicated in developmental patterning, cell-cell communication, hormonal regulation, immune function/inflammatory response, and neurotransmission. Density of DNA methylation at selected genes from the result set was also significantly associated with toddler socioemotional and behavioral problems. This is the first report to identify novel regions of the human infant genome at which DNA methylation patterns are associated longitudinally both with maternal characteristics and with offspring socioemotional and behavioral problems in toddlerhood.

Increased activation product of complement 4 protein in plasma of individuals with schizophrenia.

Structural variation in the complement 4 gene (C4) confers genetic risk for schizophrenia. The variation includes numbers of the increased C4A copy number, which predicts increased C4A mRNA expression. C4-anaphylatoxin (C4-ana) is a C4 protein fragment released upon C4 protein activation that has the potential to change the blood-brain barrier (BBB). We hypothesized that elevated plasma levels of C4-ana occur in individuals with schizophrenia (iSCZ). Blood was collected from 15 iSCZ with illness duration < 5 years and from 14 healthy controls (HC). Plasma C4-ana was measured by radioimmunoassay. Other complement activation products C3-ana, C5-ana, and terminal complement complex (TCC) were also measured. Digital-droplet PCR was used to determine C4 gene structural variation state. Recombinant C4-ana was added to primary brain endothelial cells (BEC) and permeability was measured in vitro. C4-ana concentration was elevated in plasma from iSCZ compared to HC (mean = 654 ± 16 ng/mL, 557 ± 94 respectively, p = 0.01). The patients also carried more copies of the C4AL gene and demonstrated a positive correlation between plasma C4-ana concentrations and C4A gene copy number. Furthermore, C4-ana increased the permeability of a monolayer of BEC in vitro. Our findings are consistent with a specific role for C4A protein in schizophrenia and raise the possibility that its activation product, C4-ana, increases BBB permeability. Exploratory analyses suggest the novel hypothesis that the relationship between C4-ana levels and C4A gene copy number could also be altered in iSCZ, suggesting an interaction with unknown genetic and/or environmental risk factors.

Cross-platform validation of neurotransmitter release impairments in schizophrenia patient-derived NRXN1-mutant neurons.

Heterozygous NRXN1 deletions constitute the most prevalent currently known single-gene mutation associated with schizophrenia, and additionally predispose to multiple other neurodevelopmental disorders. Engineered heterozygous NRXN1 deletions impaired neurotransmitter release in human neurons, suggesting a synaptic pathophysiological mechanism. Utilizing this observation for drug discovery, however, requires confidence in its robustness and validity. Here, we describe a multicenter effort to test the generality of this pivotal observation, using independent analyses at two laboratories of patient-derived and newly engineered human neurons with heterozygous NRXN1 deletions. Using neurons transdifferentiated from induced pluripotent stem cells that were derived from schizophrenia patients carrying heterozygous NRXN1 deletions, we observed the same synaptic impairment as in engineered NRXN1-deficient neurons. This impairment manifested as a large decrease in spontaneous synaptic events, in evoked synaptic responses, and in synaptic paired-pulse depression. Nrxn1-deficient mouse neurons generated from embryonic stem cells by the same method as human neurons did not exhibit impaired neurotransmitter release, suggesting a human-specific phenotype. Human NRXN1 deletions produced a reproducible increase in the levels of CASK, an intracellular NRXN1-binding protein, and were associated with characteristic gene-expression changes. Thus, heterozygous NRXN1 deletions robustly impair synaptic function in human neurons regardless of genetic background, enabling future drug discovery efforts.

PRG2 and AQPEP are misexpressed in fetal membranes in placenta previa and percreta†.

The obstetrical conditions placenta accreta spectrum (PAS) and placenta previa are a significant source of pregnancy-associated morbidity and mortality, yet the specific molecular and cellular underpinnings of these conditions are not known. In this study, we identified misregulated gene expression patterns in tissues from placenta previa and percreta (the most extreme form of PAS) compared with control cases. By comparing this gene set with existing placental single-cell and bulk RNA-Seq datasets, we show that the upregulated genes predominantly mark extravillous trophoblasts. We performed immunofluorescence on several candidate molecules and found that PRG2 and AQPEP protein levels are upregulated in both the fetal membranes and the placental disk in both conditions. While this increased AQPEP expression remains restricted to trophoblasts, PRG2 is mislocalized and is found throughout the fetal membranes. Using a larger patient cohort with a diverse set of gestationally aged-matched controls, we validated PRG2 as a marker for both previa and PAS and AQPEP as a marker for only previa in the fetal membranes. Our findings suggest that the extraembryonic tissues surrounding the conceptus, including both the fetal membranes and the placental disk, harbor a signature of previa and PAS that is characteristic of EVTs and that may reflect increased trophoblast invasiveness.

Comprehensive identification of somatic nucleotide variants in human brain tissue.

BACKGROUND: Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. RESULTS: Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. CONCLUSIONS: This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Glucose-limiting conditions induce an invasive population of MDA-MB-231 breast cancer cells with increased connexin 43 expression and membrane localization.

Gap junctional intercellular communication (GJIC) is a homeostatic process mediated by membrane channels composed of a protein family known as connexins. Alterations to channel activity can modulate suppression or facilitation of cancer progression. These varying roles are influenced by the cancer cell genetic profile and the context-dependent mechanisms of a dynamic extracellular environment that encompasses fluctuations to nutrient availability. To better explore the effects of altered cellular metabolism on GJIC in breast cancer, we generated a derivative of the triple-negative breast cancer cell line MDA-MB-231 optimized for growth in low-glucose. Reduced availability of glucose is commonly encountered during tumor development and leads to metabolic reprogramming in cancer cells. MDA-MB-231 low-glucose adapted cells exhibited a larger size with improved cell-cell contact and upregulation of cadherin-11. Additionally, increased protein levels of connexin 43 and greater plasma membrane localization were observed with a corresponding improvement in GJIC activity compared to the parental cell line. Since GJIC has been shown to affect cellular invasion in multiple cancer cell types, we evaluated the invasive qualities of these cells using multiple three-dimensional Matrigel growth models. Results of these experiments demonstrated a significantly more invasive phenotype. Moreover, a decrease in invasion was noted when GJIC was inhibited. Our results indicate a potential response of triple-negative breast cancer cells to reduced glucose availability that results in changes to GJIC and invasiveness. Delineation of this relationship may help elucidate mechanisms by which altered cancer cell metabolism affects GJIC and how cancer cells respond to nutrient availability in this regard.

Machine learning reveals bilateral distribution of somatic L1 insertions in human neurons and glia.

Retrotransposons can cause somatic genome variation in the human nervous system, which is hypothesized to have relevance to brain development and neuropsychiatric disease. However, the detection of individual somatic mobile element insertions presents a difficult signal-to-noise problem. Using a machine-learning method (RetroSom) and deep whole-genome sequencing, we analyzed L1 and Alu retrotransposition in sorted neurons and glia from human brains. We characterized two brain-specific L1 insertions in neurons and glia from a donor with schizophrenia. There was anatomical distribution of the L1 insertions in neurons and glia across both hemispheres, indicating retrotransposition occurred during early embryogenesis. Both insertions were within the introns of genes (CNNM2 and FRMD4A) inside genomic loci associated with neuropsychiatric disorders. Proof-of-principle experiments revealed these L1 insertions significantly reduced gene expression. These results demonstrate that RetroSom has broad applications for studies of brain development and may provide insight into the possible pathological effects of somatic retrotransposition.

Network Effects of the 15q13.3 Microdeletion on the Transcriptome and Epigenome in Human-Induced Neurons.

BACKGROUND: The 15q13.3 microdeletion is associated with several neuropsychiatric disorders, including autism and schizophrenia. Previous association and functional studies have investigated the potential role of several genes within the deletion in neuronal dysfunction, but the molecular effects of the deletion as a whole remain largely unknown. METHODS: Induced pluripotent stem cells, from 3 patients with the 15q13.3 microdeletion and 3 control subjects, were generated and converted into induced neurons. We analyzed the effects of the 15q13.3 microdeletion on genome-wide gene expression, DNA methylation, chromatin accessibility, and sensitivity to cisplatin-induced DNA damage. Furthermore, we measured gene expression changes in induced neurons with CRISPR (clustered regularly interspaced short palindromic repeats) knockouts of individual 15q13.3 microdeletion genes. RESULTS: In both induced pluripotent stem cells and induced neurons, gene copy number change within the 15q13.3 microdeletion was accompanied by significantly decreased gene expression and no compensatory changes in DNA methylation or chromatin accessibility, supporting the model that haploinsufficiency of genes within the deleted region drives the disorder. Furthermore, we observed global effects of the microdeletion on the transcriptome and epigenome, with disruptions in several neuropsychiatric disorder-associated pathways and gene families, including Wnt signaling, ribosome function, DNA binding, and clustered protocadherins. Individual gene knockouts mirrored many of the observed changes in an overlapping fashion between knockouts. CONCLUSIONS: Our multiomics analysis of the 15q13.3 microdeletion revealed downstream effects in pathways previously associated with neuropsychiatric disorders and indications of interactions between genes within the deletion. This molecular systems analysis can be applied to other chromosomal aberrations to further our etiological understanding of neuropsychiatric disorders.

Complex mosaic structural variations in human fetal brains.

Somatic mosaicism, manifesting as single nucleotide variants (SNVs), mobile element insertions, and structural changes in the DNA, is a common phenomenon in human brain cells, with potential functional consequences. Using a clonal approach, we previously detected 200-400 mosaic SNVs per cell in three human fetal brains (15-21 wk postconception). However, structural variation in the human fetal brain has not yet been investigated. Here, we discover and validate four mosaic structural variants (SVs) in the same brains and resolve their precise breakpoints. The SVs were of kilobase scale and complex, consisting of deletion(s) and rearranged genomic fragments, which sometimes originated from different chromosomes. Sequences at the breakpoints of these rearrangements had microhomologies, suggesting their origin from replication errors. One SV was found in two clones, and we timed its origin to ∼14 wk postconception. No large scale mosaic copy number variants (CNVs) were detectable in normal fetal human brains, suggesting that previously reported megabase-scale CNVs in neurons arise at later stages of development. By reanalysis of public single nuclei data from adult brain neurons, we detected an extrachromosomal circular DNA event. Our study reveals the existence of mosaic SVs in the developing human brain, likely arising from cell proliferation during mid-neurogenesis. Although relatively rare compared to SNVs and present in ∼10% of neurons, SVs in developing human brain affect a comparable number of bases in the genome (∼6200 vs. ∼4000 bp), implying that they may have similar functional consequences.

Integrated functional genomic analyses of Klinefelter and Turner syndromes reveal global network effects of altered X chromosome dosage.

In both Turner syndrome (TS) and Klinefelter syndrome (KS) copy number aberrations of the X chromosome lead to various developmental symptoms. We report a comparative analysis of TS vs. KS regarding differences at the genomic network level measured in primary samples by analyzing gene expression, DNA methylation, and chromatin conformation. X-chromosome inactivation (XCI) silences transcription from one X chromosome in female mammals, on which most genes are inactive, and some genes escape from XCI. In TS, almost all differentially expressed escape genes are down-regulated but most differentially expressed inactive genes are up-regulated. In KS, differentially expressed escape genes are up-regulated while the majority of inactive genes appear unchanged. Interestingly, 94 differentially expressed genes (DEGs) overlapped between TS and female and KS and male comparisons; and these almost uniformly display expression changes into opposite directions. DEGs on the X chromosome and the autosomes are coexpressed in both syndromes, indicating that there are molecular ripple effects of the changes in X chromosome dosage. Six potential candidate genes (RPS4X, SEPT6, NKRF, CX0rf57, NAA10, and FLNA) for KS are identified on Xq, as well as candidate central genes on Xp for TS. Only promoters of inactive genes are differentially methylated in both syndromes while escape gene promoters remain unchanged. The intrachromosomal contact map of the X chromosome in TS exhibits the structure of an active X chromosome. The discovery of shared DEGs indicates the existence of common molecular mechanisms for gene regulation in TS and KS that transmit the gene dosage changes to the transcriptome.

Epigenetic signatures of attachment insecurity and childhood adversity provide evidence for role transition in the pathogenesis of perinatal depression.

Early life adversity and insecure attachment style are known risk factors for perinatal depression. The biological pathways linking these experiences, however, have not yet been elucidated. We hypothesized that overlap in patterns of DNA methylation in association with each of these phenomena could identify genes and pathways of importance. Specifically, we wished to distinguish between allostatic-load and role-transition hypotheses of perinatal depression. We conducted a large-scale analysis of methylation patterns across 5 × 106 individual CG dinucleotides in 54 women participating in a longitudinal prospective study of perinatal depression, using clustering-based criteria for significance to control for multiple comparisons. We identified 1580 regions in which methylation density was associated with childhood adversity, 3 in which methylation density was associated with insecure attachment style, and 6 in which methylation density was associated with perinatal depression. Shorter telomeres were observed in association with childhood trauma but not with perinatal depression or attachment insecurity. A detailed analysis of methylation density in the oxytocin receptor gene revealed similar patterns of DNA methylation in association with perinatal depression and with insecure attachment style, while childhood trauma was associated with a distinct methylation pattern in this gene. Clinically, attachment style was strongly associated with depression only in pregnancy and the early postpartum, whereas the association of childhood adversity with depression was time-invariant. We concluded that the broad DNA methylation signature and reduced telomere length associated with childhood adversity could indicate increased allostatic load across multiple body systems, whereas perinatal depression and attachment insecurity may be narrower phenotypes with more limited DNA methylation signatures outside the CNS, and no apparent association with telomere length or, by extension, allostatic load. In contrast, the finding of matching DNA methylation patterns within the oxytocin receptor gene for perinatal depression and attachment insecurity is consistent with the theory that the perinatal period is a time of activation of existing attachment schemas for the purpose of structuring the mother-child relationship, and that such activation may occur in part through specific patterns of methylation of the oxytocin receptor gene.

Structural variation in the sequencing era.

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Metabolism of cholesterol and progesterone is differentially regulated in primary trophoblastic subtypes and might be disturbed in recurrent miscarriages.

During pregnancy, extravillous trophoblasts (EVTs) invade the maternal decidua and remodel the local vasculature to establish blood supply for the growing fetus. Compromised EVT function has been linked to aberrant pregnancy associated with maternal and fetal morbidity and mortality. However, metabolic features of this invasive trophoblast subtype are largely unknown. Using primary human trophoblasts isolated from first trimester placental tissues, we show that cellular cholesterol homeostasis is differentially regulated in EVTs compared with villous cytotrophoblasts. Utilizing RNA-sequencing, gene set-enrichment analysis, and functional validation, we provide evidence that EVTs display increased levels of free and esterified cholesterol. Accordingly, EVTs are characterized by increased expression of the HDL-receptor, scavenger receptor class B type I, and reduced expression of the LXR and its target genes. We further reveal that EVTs express elevated levels of hydroxy-delta-5-steroid dehydrogenase 3 beta- and steroid delta-isomerase 1 (HSD3B1) (a rate-limiting enzyme in progesterone synthesis) and are capable of secreting progesterone. Increasing cholesterol export by LXR activation reduced progesterone secretion in an ABCA1-dependent manner. Importantly, HSD3B1 expression was decreased in EVTs of idiopathic recurrent spontaneous abortions, pointing toward compromised progesterone metabolism in EVTs of early miscarriages. Here, we provide insights into the regulation of cholesterol and progesterone metabolism in trophoblastic subtypes and its putative relevance in human miscarriage.